Origins of periodic and chaotic dynamics in microfluidic loop devices 
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Droplets moving in a microfluidic loop device exhibit both periodic and chaotic behaviors based on 
the inlet droplet spacing. We propose that the periodic behavior is an outcome of a dispersed phase 
conservation principle. This conservation principle translates into a droplet spacing conservation 
equation. Additionally, we define a simple technique to identify periodicity in experimental systems 
with input scatter. Aperiodic behavior is observed in the transition regions between different periodic 
behaviors. We propose that the cause for aperiodicity is the synchronization of timing between 
the droplets entering and leaving the system. We derive an analytical expression to estimate the 
occurrence of these transition regions as a function of system parameters. We provide experimental, 
simulation and analytical results to validate the proposed theory. 



Droplets' motion in a microfluidic loop device pQ is 
a very simple phenomenon that exhibits complex spa- 
tiotemporal dynamics; such phenomena are rare in prac- 
tice [2]. A typical microfluidic loop transforms the spac- 
ing between the droplets as they travel through the de- 
vice. It has been shown that a stream of droplets en- 
tering the loop at constant inlet spacing could exit the 
loop at periodic or chaotic intervals [TJ [3HS] • The peri- 
odic and chaotic behaviors can be observed by changing 
the inlet spacing of the droplets. While the rules gov- 
erning the droplets' motion are uncomplicated, discrete 
decision-making at bifurcations and collective hydrody- 
namics lead to the intricate dynamics that are observed 
[6]. This contrast between the simplicity of the device 
and the plurality of its response has resulted in this sys- 
tem being investigated by several researchers. In all these 
studies, the broad questions of interest are: (i) what are 
the origins of the periodic and chaotic behavior? (ii) is 
there some common physics that is valid in both the pe- 
riodic and chaotic regions ? and (iii) is it possible to 
predict the regions of periodic and chaotic behavior? Al- 
though the existing literature partially answers some of 
these questions such as the prediction of inlet spacings 
where transitions from periodic to chaotic behavior oc- 
cur HI |5j [3 E] , the cause for these transitions is not well 
understood. In this letter, we seek to address the ques- 
tions posed above as comprehensively as possible at this 
time. While understanding the loop dynamics is in it- 
self intellectually compelling, the insights gleaned from 
such an exercise could also help in understanding natu- 
ral systems such as blood flow in microvasculature and 
in developing practical design and control approaches in 
the field of discrete droplet microfluidics. 

Consider a train of droplets entering a microfluidic loop 
device as shown in Fig [I] at a constant inlet spacing A. 
The relative spacings between the droplets remain con- 
stant until they reach the loop junction. As the droplets 
enter the upper or lower branch, the resistance of that 
branch changes. This results in a redistribution of the 
bulk flows to the upper and lower branches. Since the 
droplets' velocities are proportional to the bulk flow, the 




FIG. 1. A constant stream of droplets passing through a 
microfluidic loop device 



droplets in the upper and lower branches move at differ- 
ent velocities leading to changes in the relative distances 
between them as they exit the loop. Mathematically, 
the loop device (G) transforms the inlet spacing (A) to 
a sequence, {A} n = {G(A)} n . It has been shown in the 
literature that the exit spacings can either form a peri- 
odic or a chaotic pattern. These patterns are captured 
through a Poincare map, where the number of distinct 
clusters represent the number of periods. There are two 
disadvantages to such a characterization: it is difficult to 
exactly identify the number of periods when confronted 
with scatter in the data, and the plot does not explain 
the physics of the process that leads to these periodic 
or aperiodic behaviors. To address these two issues, we 
start with a mathematical definition of periodicity. The 
exit sequence ({A} n ) is periodic with period p, where p 
is the smallest integer such that = {A}i +P V i e N. 
The exit sequence is aperiodic, if there exists no p such 
that = {\}i+ p V i e N. From the definition of pe- 
riodicity it can be noticed that the periodicity condition 
will hold for all integer multiples of p. Further, it is self 
evident that the periodicity condition cannot hold for 
two p's that are not multiples of each other. However, 
if the periodicity condition is valid in some region, i.e., 
i e (iVi,iV2) for some p, then that behavior is termed 
locally periodic. We are interested in the physics that 
leads to periodic behavior and the manner in which the 
system transitions from periodic to chaotic behavior and 
vice versa. 

The sustained periodic behavior is a result of a series 
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of droplets taking the same decisions repeatedly - this 
can also be viewed as memory of the system [8 j. The 
decisions will be repeated with period p if the i th and 
i + p th droplets undergo the same experience i.e., veloc- 
ity changes as they pass through the loop. This hap- 
pens if the loop snapshot (by which we mean the device 
with droplets at particular positions in the loop) that 
the i + p th droplet encounters is the same as the one en- 
countered by the i th droplet. A natural consequence of 
this concept is that the amount of the bulk phase fluid 
trapped between the i th and i + p th droplets should ex- 
actly equal the amount of the fluid trapped between the 
same droplets prior to entering the loop. Mathematically, 
the bulk fluid present in between i th and i + p th droplets 
prior to the device entry is pXS, where S is the cross 
sectional area of the channel. Assuming bulk phase in- 
compressibility, the volume trapped between the droplets 
after exit is the sum of all the exit relative spacings mul- 
tiplied by S. This volume conservation that is required 
for repetition translates to a remarkable "conservation of 
droplet spacing" given below: 
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If there exists a unique p, V i e N then eqn[T]is an 'if and 
only if condition for a loop device to show sustained pe- 
riodic oscillations. We validate this 'conservation princi- 
ple' that is derived through physical arguments using the 
network model proposed in pQ . We model two experimen- 
tal scenarios, 3-period and 9-period behaviors, without 
input fluctuations. We plot the value of |S^ _1 A— pX\/X 
as a function of j as shown in Fig [2] If our insight is cor- 
rect then this value should be equal to zero whenever j is 
equal to multiples of p, the periodicity of the system. The 
initial (I) that is used in the equation does not have any 
impact on the result because the case studies considered 
here have sustained periodic oscillations. As expected, in 
Fig [2j the curve touches the x-axis at N = 3,6,9... for 
the 3-period case and N = 9, 18 ... for the 9-period case, 
validating the 'conservation principle'. 




FIG. 2. Conservation of droplet spacing in the case of 3 and 
9 period behaviors 

Estimation of the period p in case of experimental 



systems can be complicated due to fluctuations in the 
inlet spacing. In such cases using conventional methods 
such as Poincare maps [6 to obtain periodicity may be 
misleading. The conservation principle provides another 
view of periodicity. Looking at eqn [I] it is apparent 
that the periodicity is equal to the number of droplets 
that exit before a snapshot repeats. Consider a 2-period 
behavior as identified by a Poincare map depicted 
in Fig [3] Modeling the system without experimental 
fluctuations at the inlet using a network model leads 
to a 3-period behavior. To confirm if this is the true 
period of the system we used the snapshot comparison 
technique derived from the conservation principle. 
The 1 st image is repeated after the 132 nd and 264 th 
snapshots as shown in the inset of Fig [3j with the 
droplets exiting at the 1 st , 132 nd and 264 th images. In 
this case, three droplets exit before the original image 
is repeated, confirming that this indeed is 3-period 
behavior, contrary to what Poincare map suggests. 
The proposed snapshot-comparison method is therefore 
particularly helpful in estimating the periodicity of loop 
dynamics where experimental scatter might otherwise 
make such identification problematic. 



1.5 



A 

I— 

v 



0.5- 



+ Simulation 
■ Experiment 



Mtfi' p .... 





132 nd 

1 




264 th 




0.4 0.6 0.8 1 1.2 1.4 1.6 

T n /<T> 



FIG. 3. Three droplets exit the loop device in between the 
images I s *, 132 nd and 264 th . Poincare map of the exit times 
shows a two period behavior, whereas the simulation shows 
three periods. Oil and water flow rates are 550 and 25 p L/hr 
respectively 

We now turn our attention to the understanding of 
aperiodic behavior and the physics that leads to the tran- 
sitions from periodic to aperiodic behavior and vice versa. 
We answer this transition question by exploring what 
happens to the conservation principle in the aperiodic 
region. There are two possibilities: one in which the 
conservation principle is never valid for any two pairs of 
droplets and the other scenario is where the conservation 
principle is valid between some droplet pairs but not re- 
peatedly ad infinitum. To resolve between the two possi- 
bilities we again take recourse to the the network model 
PQ where we simulate a microfluidic loop for aperiodic 
behavior. In Fig [4j the x-axis of the plot is increasing 
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natural numbers (/) representing the I th time that the 
conservation principle holds and the y-axis representing 
the number of the droplets that participate in the con- 
servation equation (the measure of periodicity) when it 
holds for the I th time. This plot suggests that even in the 
aperiodic case, the conservation principle is still driving 
the physics; however, the periodicity changes with time 
and the conservation equation is valid at all these peri- 
odicities. Inset Fig [4^ a) shows the initial period, where 
for a while a local periodic behavior is observed, which 
is broken, followed by the same period for a shorter time 
and so on. This is a classic case of intermittent behavior; 
similar behavior is observed in inset Fig[4^b). Inset Fig 
[4^c) shows the dynamics settling down to a three period 
behavior for a larger duration but the conservation prin- 
ciple holds for alternate 1 and 2 droplets and, as a result 
for 3 droplets at all times, demonstrating the richness of 
the behavior exhibited by this simple loop device. Fig 
[4^d) shows an experimental system displaying aperiodic 
behavior where the conservation principle is valid for dis- 
tinct sets of droplets. This result has implications on how 
the system transitions from periodic to aperiodic behav- 
iors and suggests that the aperiodic behavior is likely to 
be interspersed between the periodic regimes. 

The intermittent transitions between periodic behav- 
iors in an aperiodic scenario is depicted in Fig [5] where 
the values '0' and '1' represent the binary decisions of 
the droplets i.e., choice of either the upper or the lower 
branch. In the transition regions any combination of 0's 
& l's is possible. This can result in large intermittent 
periods and sporadic transitions to other periodic (for 
example 2 or 3) behaviors. It can be noticed that a re- 
markably similar picture is observed in the intermittency 
graph derived from the conservation principle as shown in 
Fig [4] Interestingly, similar phenomena are also observed 
in cellular automata [9]. Fig|6]shows an experimental sce- 
nario where the intermittent behavior is observed. Con- 
sider Fig|6ja); the snapshots are arranged in an increas- 
ing order of time. In all these images the droplets enter 
the loop when the loop contains 5 droplets. As the exiting 
drop leaves, the entering drop can either take the lower 
or the upper branch. In the case of Fig[6^a), the enter- 
ing droplet chooses the upper branch, but in the case of 
Fig[6jb), the droplet chooses the lower branch. Although 
the system has very similar conditions at both times, the 
dynamics of this system moves to a completely different 
regime due to a single droplet altering its decision. Since 
such transitions happen at different time instants during 
the course of the droplets' travel through the loop, the 
dynamics of the system becomes unpredictable. 

In summary, the transition region is where the con- 
servation principle holds for several periods at different 
times; if the droplets take decisions in a repetitive fashion 
ad infinitum, then the behavior is periodic. For aperi- 
odic behavior to be manifested, there has to be a critical 
event to break the constant periodic cycle. We propose 



that this critical event is one where the droplets enter and 
exit at around the same time, leading to changes in the 
droplets' decisions. This will in turn change the amount 
of dispersed phase fluid trapped between the droplets. 
Predicting the onset of transitions is of significant impor- 
tance. Reliable prediction of the transition region can be 
used to design experiments and to manufacture devices 
that will consistently show either periodic or chaotic be- 
havior. It is clear that there are an infinite set of A val- 
ues in probably a very small range that characterizes the 
transition region. Predicting the range of A values where 
the transitions occur is a difficult problem. Instead, in 
this work, we attempt to predict a particular value of 
A (A c ) in the transition region. From a practical stand- 
point, since we expect the transition regions to be char- 
acterized by a very small range of A values, prediction 
of one value in the range should suffice. To predict one 
A c , we need to find a A that results in a scenario where 
a droplet enters the loop device while another droplet 
is exiting. While this could happen for several different 
pairs of droplets (leading to a range of A values), we as- 
sume that the critical A c is when the first droplet exits 
the loop device while another droplet enters the device. 
However, depending on the A c , the number of droplets 
that a loop holds (n) before the first droplet exits could 
be different. Remarkably, we will demonstrate that by 
changing n, all the transition regions can be identified. 
Interestingly, the n used in the calculation of a particu- 
lar transition A c does not provide any indication of the 
period p of the sustained periodic behavior that follows 
the transition. This is because the transition A c has no 
bearing on the number of droplets for which the con- 
servation principle holds. Consider Fig{l]- the entering 
droplets move at a velocity v = /3Q/S, where j3 is the 
slip factor, Q, S, L\ and L u are inlet flow rate, cross sec- 
tional area of the channel and length of upper and lower 
branches respectively. Assuming that the loop is sym- 
metric (Li & L u ), droplets alternate between top and 
bottom branches. The resistance of the lower branch is 
R = Ri + mRd, where Ri is the branch resistance, Rj is 
the droplet resistance and m is the number of droplets in 
the branch. 

The exit time of the first droplet after n droplets enter 
the loop device is as follows: 
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In the above equation, t\, ex n is the time taken by the 
droplet from that instant to exit the loop, as shown in 
Fig [I] and is not the total time it takes to exit the loop. 
The entry time of the next droplet in the loop is given 
by t entry = ^| • Therefore the critical (A c ) values occur 
when: t\^xit — Gentry This results in (with R = R^/ R u ) 



4 



70 
60 
50 

I 40 

O 

6 30 
z 

20 
10 
Q 




200 




(b) 


2 






100 






1.5 
1 


h/w\ 


/wv 








5 
















40 50 60 290 300 310 305 310 315 



m (a) l 



h 








iTff 






150^ 200 



FIG. 4. Periodicity calculated using Eqn[T]for an aperiodic 
system at distinct intervals of time index; (a),(b),(c) show 
distinct regions of the above plot, inset (d) shows an aperiodic 
experimental system demonstrating conservation of spacing. 
Oil and water flow rates are 250 and 25 fi L/hr respectively 
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FIG. 5. Representation of transition phenomena in a microflu- 
idic loop device 




FIG. 6. A loop device showing transitions: (a) A sequence 
of snapshots showing the entering droplet choosing the upper 
branch; (b) At a later time in the same experiment, with sim- 
ilar conditions, the entering droplet chooses the lower branch. 
Oil and water flow rate are 250 and 25 ji L/hr respectively 
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The critical A c values are obtained by varying n, n G N. 
The proposed analytical expression is validated using a 
bifurcation map along with an analytical expression from 
the literature [4|. 

A particular loop configuration is simulated to validate 
the theoretical predictions as shown in Fig [7| The flow 
rate is kept constant, and only the feeding frequency is 
varied. Based on input frequency, after the initial tran- 
sients, the droplet exit pattern reaches a steady state. At 
this point, periodicity is calculated and plotted against 
the inlet spacing. Fig. [7^a) shows a bifurcation map, 
where normalized inlet spacing is plotted on the X-axis 
and the normalized outlet spacing is plotted on the y- 
axis. 

At A > 1, the system shows a steady 2-period behav- 
ior, but as the spacing decreases, the outlet periodicity 
transitions from a 2-period to a 3-period behavior - ob- 
serve that there is no period doubling vis a vis logistic 
map [TO] . The 3-period behavior lasts for some time and 
the system goes back to a two period behavior. After 



FIG. 7. Bifurcation map simulated by changing the input 
spacing in a scenario where the oil and water flow rates are 
550 and 25 \i L/hr respectively, the bottom plot shows the 
prediction of transitions in the bifurcation maps using the 
proposed analytical expression along with that reported in 
the literature [3]. The upper and lower branch lengths of the 
loop device are 2235/xra and 2217/ira, respectively 

the next transition, the system drifts to a 4-period be- 
havior. We observe only 3 values repeating within the 
4-period pattern. The analytical expression captures all 
these transition regions corroborating the intuition that 
droplets entering and exiting at the same time leads to 
chaotic behavior. The results from the expression for 
critical A c values proposed in the literature are also plot- 
ted in Fig. [7^b) and are in agreement with the theory 
that is proposed in this paper. 

In conclusion, in this paper, we propose a conservation 
principle that is valid for periodic and chaotic behav- 
ior in microfluidic loop systems. Our experimental and 
theoretical analysis provides a unique way to quantify 
this behavior. As the conservation of spacing is valid for 
both periodic and chaotic regions of the loop device, this 
device resembles a Hamiltonian system [11 . The chaotic 
behavior is a result of transitions occurring when droplets 
enter and exit at the same time. This is akin to complex 
sequences derived with simple rules in Cellular Automata 



[9] . We derive an analytical expression for predicting the 
transition regions based on the new physical insights pro- 
posed in this paper. The results suggest that the system 
achieves chaotic behavior through a phenomenon called 
"intermittency route to chaos" [lOj [12] and not just be- 
cause of a 3 period behavior [T3] . 
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